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Abstract. We report the discovery of a Neptune-mass planet around Gl 581 (M3V, M = 0.31 M G ), based on precise Doppler 
measurements with the HARPS spectrograph at La Silla Observatory. The radial velocities reveal a circular orbit of pe- 
riod P = 5.366 days and semi-amplitude Ki = 13.2 m s -1 . The resulting minimum mass of the planet (w^sini) is only 
0.052 Mj up = 0.97 M Nep = 16.6 M E „ rth making Gl 581b one of the lightest extra-solar planet known to date. The Gl 581 
planetary system is only the third centered on an M dwarf, joining the Gl 876 three-planet system and the lone planet around 
Gl 436. Its discovery reinforces the emerging tendency of such planets to be of low mass, and found at short orbital periods. 
The statistical properties of the planets orbiting M dwarfs do not seem to match a simple mass scaling of their counterparts 
around solar- type stars. 
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1. Introduction 

Over 150 planets have been found orbiting main se- 
quence stars other than the Sun, in about 140 plan- 
etary systems of which 18 have multiple planets 
(http://vo.obspm.fr/exoplanetes/encyclo/!. These extra- 
solar planets are a very diverse class: their mass ranges 
between half the mass of Neptune and 15 times the mass 
of Jupiter, some have large eccentricities when others have 
nearly circular orbits , their periods range from slightly over 
a dav jK onacki et "ail 120031 OGLE-TR-56) to over a decade 
(Mar cv et all l2002t 55 Cnc). The multiple s ystems range 
from strongly resonant to fully hierarchical jRivera et all 
l2005UMarcvetali r2002). This diversity demonstrates that our 
own solar system represents but one possible outcome of the 
planetary formation and evolution processes, and apparently 
not even a very common one. 

The statistical properties of these exoplanets provide cru- 
cial clues to their formation mechanism. As perhaps the most 
dramatic exampl e, the seminal detection of the 51 Peg planet 
in a 4-days orbit IMavor & Ouelozul9 95) immediately forced 
theoreticians to recognize the critical importance of orbital mi- 



* Based on observations made with the HARPS instrument on the 
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gration dLin et allll99ritlWardLfl997l) . The correlation between 
the occurence of Jupiter- mass planet and the high metall icity 
of the host stars JGonzaleztri997HSantos et all 1200 lll2004 is 
another exemple. It is thought to reflect the controlling role of 
the condensate mass in the protoplanetary disk, but it has taken 
longer to converge towards that consensus. 

To date, all but 2 of these 140 planetary systems or- 
bit solar-type stars. In part, this no doubts reflects a bias of 
most planet-search programmes towards the relatively bright 
F to K main sequence stars, and away from their fainter M- 
type counterparts (M < 0.6 Mb ). Non etheless, severa l teams 
dWright et all Eool iRndl et all Eool lOelfosse e"Hil M 998ah 
collectively monitor over 200 M dwarfs with sufficient pre- 
cision to detect a Jupiter-mass planet out to at least 2 AU. 
These efforts have up to now identified the 3-planet sys- 
tem around Gl 876 dDelfosse et all Il998al iMarcv et all Il998l 
2001; Rivera et al., 2005), and the single-planet Gl 436 system 
(But ler et all 12004). Of these 4 planets, 2 are in the Neptune- 
mass class, leaving only two of the Gl 876 planets with ap- 
proximately Jupiter-mass. By constrast > 5% of solar-type stars 
have Jupiter-mass planets ( Mar cv etall l2000l) . and the com- 
para tive deficit for the M dwarfs is the refore statistically robust 
dButler et alll2004lNaef et al. L 1200.4 . 

An open question, though, is whether M dwarfs gen- 
uinely have fewer planets, or whether their planets are just 
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as abundant, but not quite as massive. Addressing this ques- 
tion needs higher precision, or more measurements, than the 
radial-velocity surveys have achieved to date. To help answer- 
ing this question, we are using the HARPS spectrograph for a 
high-precision survey of more than 100 nearby M dwarfs. We 
present here its first detection, a Neptune-mass planet around 
G1581. 

2. Properties of Gl 581 

G1581 (HIP 74995, LHS 394) is an M3 dwarf JHawlev et all 



1997|) with a distance to the Sun of 6.3 pc (n = 159.52+2.27, 
lESAl Jl997ll). Its phot ometry (V=10. 55+0.01, B-V =1.60 
llMermimodetailll997l) : K=5.85±0.03 dLeggettl fl99l l and 
the parallax together result in absolute magnitudes of 
My=11.56±0.03 and M K =6.86+0.04. From its absolute V 
magnitu de and the 2.08 V-band bolometric correction of 
belfosse et all lll998h. the luminosity of Gl 581 is 0.013 L . 
The iDelfosseetalJ (I2000h K-band mass-luminosity relation, 
which has much lower intrinsic dispersion than the equiv- 
alent V-band relation, giv es a 0.31+0.02 M mass, and the 
Chabrier & Baraffe] fcOOOh theoretic al Mass-Radius rel ation 
then a radius of 0.29 R . Interestinglv. lBonfils et alj f2005 ) find 
Gl 581 slightly metal-poor ([Fe/H]=-0.25), in contrast to most 
planet-host stars having supersolar metallicities. 

Gl 581 has been classified as a variable star (HO Lib) . 
However, the data which have led to this classification (Weis, 
1994) have a short-term variability of ~ 0.006 mag. The vari- 
ability quoted by the author is marginally above the errors bars 
and, if real, has most likely a long-term nature (several years). 

The age of Gl 58 1 can be estimated from its kinematic char- 
acteristics, its magnetic activity, and its metallicity, all of w hich 
point towards a moderate to older age. Leggett ( 1992) find 
that its UVW galactic velocities are intermediate between those 
typica l of the young and old galactic disk, and Del fosse et alJ 
Jl998h find very low X-ray emission (L x /Lt, i < 5.10~ 6 ) and a 
2. 1 km s _1 upper limit on the projected rotation velocity, v sini. 
The HARPS spectra show weak Ca// H and K emission, in 
the lower quartile of stars with similar spectral types (Fig.Q. 
As mentioned above, Gl 581 also has a subsolar metallicity. 
Altogether, these properties suggest that it is at least 2 Gyr old, 
and they ensure that the radial velocity "jitter" from magnetic 
activity must be minimal. 

3. Doppler measurements and orbital analysis 

HARPS (High Accuracy Radial velocity Planet Searcher) is 
the new ESO high-resolution (R = 115 000) fiber-fed echelle 
spectrograph, optimised for planet search programmes and as- 
teroseismology. It has proved to be the most precise spectro- 
velocimeter to date, reaching an instrumental RV accuracy 
better than 1 m s" 1 JMavor et all 120031 ISantos et all |2004a; 
IPepe et all 120041 Lovis et all 120051) . and even better on the 
short-term scales of interest for asteroseismology. For ultimate 
radial velocity precision HARPS uses simultaneous exposures 
of a thorium lamp through a calibration fiber. When observ- 
ing M dwarfs however, we rely instead on its very high instru- 
mental stability (nightly instrumental drifts < 1 m s~'). The M 
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Fig. 1. HARPS spectra of the Call H (A = 3968.47 A) line re- 
gion for Gl 58 1 and three comparison stars with similar spectral 
type and apparent magnitude. The stars are displayed in order 
of ascending chromospheric activity, and from top to bottom 
are G1752A, G1908, G1832 and G1581. The chromospheric 
emission peaks look prominent against the weak blue contin- 
uum of these M dwarfs, but they actually denote very weak 
chromospheric emission relative to the bolometric luminosity. 
Amongst those 4 stars, G1581 has the weakest chromospheric 
activity. 



dwarfs are typically too faint for us to reach the stability limit of 
HARPS within realistic integration times, and dispensing with 
the simultaneous thorium light produces much cleaner stellar 
spectra, suitable for quantitative spectroscopic analyses. 

For the V = 10.5 Gl 581 we use 15 mn exposures, and 
the median S/N ratio of our 20 spectra is 40 per pixel at 
550 nm. The radial velocities (Table |2] only available elec- 
tronically) were obtained with the standard HARPS reduction 
pipeline, based on the cross-correlation with a stellar template 
and the precise nightly w avelength calibration with ThAr spec- 
tra JBaranne et allfl996l) . They have a median internal error of 
only 1.3 m/s, which includes both the nightly zero-point cali- 
bration uncertainty (~ 0.8 m s~') and the photon noise, com- 
puted fro m the full D oppler information content of the spectra 
dBouchv et"aill200lh . 
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Fig. 2. Upper panel: Phased radial velocities for Gl 58 1 . Lower 
panel: Residuals around the fitted solution versus time. 
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Fig. 3. Bissectors of the HARPS correlation profiles for Gl 58 1 . 
The shape of the bissector curve does not change together with 
its position. It validates the planetary interpretation as the cause 
of the radial-velocity variation. 

Table 1. Orbital and physical parameters. 



The computed velocities exhibit an rms dispersion of 
10 m s _1 , much above their internal errors and also consider- 
ably more than we observe for stars with higher chromospheric 
activity. Of the three comparison stars with stronger chromo- 
spheric emission in Fig. [2 G1752 and G1832 have enough 
HARPS radial velocities to measure rms dispersions of 1.2 m/s 
(from 10 measurem ents) and 1.8 m/s (from 19 measurements). 
iRivera et al.l J2005I) report an rms dispersion of 3.42 m/s for the 
the third, G1908, dominated by their measurement noise. 

As demonstrated by Fig|2] a circular orbit of period 
P = 5.366 d and semi-amplitude K = 13.2 m s _1 is an excel- 
lent fit to these velocities, as expected from rapid tidal cir- 
cularisation at this short period. Attempts at adjusting ellipti- 
cal orbits resulted in non-significant eccentricities. We there- 
fore adopt a circular orbit. Its parameters are listed in Tabled 
and lead to a minimum mass (M sin i) for the planet of only 
0.052 M Jup = 0.97 M Nep = 16.6 M Earlh . The weighted rms of 
the residuals around the fit is 2.5 m s _1 , and twice the internal 
errors of the measurements. More data points are needed to es- 
tablish whether this extra dispersion is intrinsic to the star, or 
whether it could be explained by the presence of a third body 
in the system. Further credit is given to the latter hypothesis 
by the very-low variability level in the shape of the bisector 
(Fig |3- An overall translation dominates the evolution of the 
spectral profile, leaving no doubt on the keplerian origin of the 
radial-velocity variations. We also have available 1 1 ELODIE 
spectra, which extend the measurement baseline to 9 years, al- 
beit with a 5 years gap between mid-2000 and mid-2005. Their 
17 m/s median error bars are too large to reveal the K=13 m s _1 
planetary signal. They exhibit an rms dispersion of 23 m s _1 
and lack any obvious long term trend, demonstrating that the 



Parameter 


G1581 b 


P 


[days] 


5.366 ± 0.001 


T 


[JD-2400000] 


51004.30 ±0.06 


e 




0. (fixed) 


V 


[kms- 1 ] 


-9.212 ± 0.001 


CO 


[deg] 


0. (fixed) 


K 


[m s- 1 ] 


13.2 ±0.4 


a\ sin i 


[10- 3 AU] 


0.652 10~ 2 


f(m) 


[10" 9 M G ] 


0.1354 10~ 2 


iri2 sin i 


[^Earlh] 


16.6 


a 


[AU] 


0.041 


N 




20 


5 pan 


[days] 


440 


a (O-C) 


[m s~'] 


2.5 



Gl 581 system contains no Jupiter-mass planet with any period 
shorter than ^10 years. 

4. Discussion 

The semi-major axis of the planetary orbit is only 0.042 AU or 
9 solar radii, similar to most close-in planets around solar-type 
stars. In the natural length unit of its central star however, this 
amounts to 3 1 radii of Gl 58 1 . The geometric transit probability 
is thus only 3%, and significantly less than the »10% typical 
of close-in planets around solar-type stars. If transits do oc- 
cur on the other hand, the planet will cover a larger fraction of 
its smaller star. For a constant planetary radius, transits would 
thus be correspondingly deeper and more easily detected. At 
this radius and given the 0.013 L luminosity of the star, the 
expected temperature of the planetary surface is ^420 K, with 
large uncertainties from the unknown albedo and energy trans- 
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port. Even with conservative error bars though, this tempera- 
ture is compatible with either a rocky planet or a gas giant, and 
evaporation will be negligibly small in either configuration. 

The detection of Gl 581b brings the inventory of M-dwarf 
which harbor planetary systems to 3, and the number of their 
planets to 5. While admittedly still very small, these samples 
allow us an initial peek at their properties as a population, 
compared to the much more numerous planets known around 
solar-type stars. One immediate observation is that none of 
the three stars is metal-rich, with Gl 876, Gl 436 and Gl 581 
having metall icities of respectiv ely [Fe/H]=-0.03, +0.02 and 
-0.25 dex JBonfils et all l2005h . This contrasts with the me- 
dian metallici ty for solar- t ype s tars surrounded by planets, 
[Fe/H]=+0.2 (Santos et al., 2005), though the significance of 
the difference is obviously still modest. 

As discussed in the introduction, various groups monitor 
over 200 M dwarfs with 3-15 m s _1 precision, sufficient to 
easily detect the >40 m s reflex motion of a 0.3M Q star 
orbited by a Jupiter-mass planet out to 2 AU. That these ef- 
forts have to date found only 5 planets, of which only Gl 876b 
and c have approximately jovian masses, demonstrates that 
there are much fewer «M.r„„ planets around M dwarfs than 
the *5% JMarcv et all l200ol iNaef et al. L 12005) found around 
solar-type stars. The 5 planets include no hot- Jupiter, but with 
only ~1% solar-type star orbited by such a planet the sig- 
nificance of that fact is still modest. 3 of the 5 on the other 
hand are hot-Neptunes (Gl 436b, M sin i=1.2 M Nep ; Gl 876d, 
0.44 M Nep ; Gl 581b, 0.99 M Nep ), as many as currently known 
arou nd all so l ar-typ e stars. This matches the theoretical model 
of llda & Linl J2005I) : the mass-distribution of close-in planets 
has two peaks centered at about the masses of Jupiter and 
Neptune, with the former preferentially populated around G- 
dwarfs and the latter around M dwarfs, reflecting how much 
matter remains available in the disk for accretion during the 
inward migration of the planet. Other theoreticians however 
take the vie w that many h ot-Neptunes are actually evaporated 
hot-Jupiters (Bara ffe et all l2005). Better statistics on M-dwarf 
planets will help determine which of these mechanisms domi- 
nate. 

A final striking characteristic of the current M-dwarf plan- 
ets is that none has a period longer than the 2 months of 
Gl 876b. By contrast, 66% of the 164 planets known around 
solar-type stars have orbital periods above 100 days and 
their distribution is even observed to increase with period 
ftldrv et all 120031) . With 5 planets known around M dwarfs, 
the probability that the long-period deficit amongst M-dwarf 
occurs by chance is thus less than (0.34) 5 =5.10~ 3 . The sensi- 
tivity of Doppler searches does degrade for longer periods how- 
ever, and planets of M dwarfs have often been found close to 
the sensitivity floor of their respective discovery surveys, mak- 
ing detection biases potentially important. A full account is not 
currently possible from published information (our own high 
precision M-dwarf survey is still too recent to be very useful in 
this respect), but for periods shorter than the observing interval 
the sensitivity degrades only slowly (ocP-1/3). Some of the M- 
dwarf Doppler surveys have been obser ving for long enough 
jBonfils et all l2004 iButler et all 2004) that detection biases 
make an unlikely full explanation for the lack of long period 



planets. That deficit therefore has to be at least in part intrin- 
sic, perhaps reflecting smaller protoplanetary disks around the 
lower mass M dwarfs. 
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Table 2. Radial-velocity measurements and error bars for 
G1581. All values are relative to the solar system barycenter. 
Only available electronically 
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